The Crab Nebula as a standard candle in 
very high-energy astrophysics 
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The continuum high-energy gamma-ray emission between 1 GeV and 10 3 GeV from the Crab Nebula has 
been measured for the first time in overlapping energy bands by the Fermi large-area telescope (Fermi/LAT) 
below ss 100 GeV and by ground-based imaging air Cherenkov telescopes (IACTs) above ~ 60 GeV. To fol- 
low up on the phenomenological approach suggested by Hillas et al. (1998), the broad band spectral and 
spatial measurement (from radio to low-energy gamma-rays < 1 GeV) is used to extract the shape of the elec- 
tron spectrum. While this model per construction provides an excellent description of the data at energies 
< 1 GeV, the predicted inverse Compton component matches the combined Fermi/LAT and IACT measure- 
ments remarkably well after including all relevant seed photon fields and fitting the average magnetic field to 
B = (l24 ± 6 (stat.) ^ 5 (sys.)) pG. The close match of the resulting broad band inverse Compton component 
with the combined Fermi/LAT and IACTs data is used to derive instrument specific energy-calibration factors. 
These factors can be used to combine data from Fermi/LAT and IACTs without suffering from systematic un- 
certainties on the common energy scale. As a first application of the cross calibration, we derive an upper 
limit to the diffuse gamma-ray emission between 250 GeV and 1 TeV based upon the combined measurements 
of Fermi/LAT and the H.E.S.S. ground-based Cherenkov telescopes. Finally, the predictions of the magneto- 
hydrodynamic flow model of Kennel & Coroniti (1984) are compared to the measured SED. 

Keywords: ISM: individual objects: Crab Nebula - Radiation mechanisms: non-thermal - Acceleration of particles - Magne- 
tohydrodynamics (MHD) - Gamma rays: diffuse Background 



I. INTRODUCTION 

The Crab Nebula has been and remains an intensel y stud- 
ied ob ject in astrophysics (for a recent review see e.g. lHesterl 
120081) . It is part of the remnant of a core-collapse super- 
nova that occu rred in 1054 AD at a distance of d ~ 2kpc 
dTrimblelll968h . Observations of the nebula have been car- 
ried out at every accessible wavelength resulting in a re- 
markably well-determined spectral energy distribution (SED). 
Therefore, the Crab Nebula is an ideal object for detailed 
studies of th e conversion of Poynting flux to particle energy 
flux (see e.g. lAronsll2008l:[Coronitilll99filKirk & SkiaeraaserJ 
12003 ) and finally of the acceleration processes taking place 
at th e termination shock (see e .g. lEmmering & Chevalier! 
119871: ISpitkovskv & Aronsl |2004|) . The com monly consid- 
ered model for the Crab Nebula (see e.g. iRees & Gunnl 
1 1974 henceforth RG74) assumes an ultra-relativistic out- 
flow from the pulsar that terminates in a standing shock at 
a distance r s which is roughly 10% of the total nebula's 
size. In the downstream medium, the particles are pitch- 
angle isotropized, forming a broad power law in energy. The 
magneto -hydrodynamic (MHD) an alysis of the downstream 
flow by iKennel & Coronitil dl984l) . henceforth KC84, pro- 
vides an elegant solution to the particle distribution and the 
magnetic field downstream of the shock under the assumption 
of a particular ad hoc injection spectrum. In the framework 
of this model, the radio emission is explained by a separate 
electron population, which has been link ed to the high spin- 
down phase of the pulsar (relic ele ctrons, Atovanll 19991) or by 
acceleration in MHD turbulences dNodes et alj|2004 . 

Several different models have been proposed to explain 
the observed high-energy gamma-rays as inverse Compton 



emission from the same electron populatio n responsible for 
the synchrotron X-ray emission (see, e.g.. lAharonian et al] 
2004-lAtovan & Aharonianll996l:lBednarek & Bartosik|2003 : 
de Jager & Hardindl992l:lHillas et alJl998tlVolpi et alj2008 : 



Zhang et all 120081 We point out that IVolpi et al] use a time 
dependent axisymmetric numerical simulation of the nebula's 
evolution). The large uncertainties in the observational data 
in the past have only weakly constrained the models at high 
energies, i.e., between 1 GeV and 100 GeV. The largely im- 
proved statistics of the measureme nts carried out with t he re- 
cently commissioned Fermi/LAT dAtwood et al] l20Q9h have 
provided us with more accurate data in this crucial energy 
window. Nevertheless, the combined observations carried 
out with different instruments in overlapping energy bands 
(Fermi/LAT and IACTs) are currently limited by the system- 
atic uncertainties of the relative and absolute energy calibra- 
tion. 

Because of this, a two-pronged approach is followed here. 
On the one hand, we perform accurate modeling of the avail- 
able measurements of the broad band SED of the Crab Neb- 
ula. On the other, the model is used to derive corrections 
of the measured energy scale for the individual instruments 
to a common energy scale. This cross calibration reduces 
the systematic uncertainties and proves useful for any study 
which relies on the combination of spectral measurements of 
the Fermi/LAT and IACTs. 

Given that we are mainly interested in the spectral modeling 
of the high-energy emission, which is currently not spatially 
resolvable, we primaril y apply a simple a nd robust approach 
based upon the work of lHillas et al. I dl998l) . which we refer to 
as the constant B-field model. Nevertheless, we also investi- 
gate the MHD flow model suggested by KC84 and its applica- 
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tion to the inverse Compton component jAtovan & Aharonianl 
1 19961 henceforth AA96). 

The article is organized as follows. In Section[n]the avail- 
able data and modeling of the SED are presented and dis- 
cussed. Both, the simplified constant B-field model and the 
MHD flow model of KC84 are fitted to the measured syn- 
chrotron part of the SED by means of a ^-minimization. The 
results of the two models are compared. In Section [III] we 
discuss the properties of the underlying electron spectrum of 
the simplified model as derived from the observations. The 
results of a first instrumental cross calibration are presented 
in Section[IV] together with an application to extract limits on 
the diffuse y-ray background at TeV energies. 



TABLE I. Observations used for th e SED. All other data s hown 
in Fig. |l(a)|(b)| are ta ken from | Aharonian et all d2004|) and 
references therein. (1 ) iGreenet alT d2004); (2) iTemim et all 



j2006h: (3)lKirschet al.N2005l): (4)|jourdain & Roquesl d2009h: (5) 



iJourdain et alj d2008h: (6) lAbdo et alj | |201C|) ; (7) lAharonian et all 



d2006l) : (8) lAlbert et alj j2008h 



Energy Band Instrument 


Reference 


Submillimeter ISO & SCUBA 


(1) 


to far infrared 


SPITZER 


(2) 


X-ray to 


XMM-Newton 


(3) 


y-ray 


SPI 


(4) 




IBIS/ISGRI 


(5) 




Fermi / LAT 


(6) 


VHE 


H.E.S.S. 


(7) 




MAGIC 


(8) 



II. SPECTRAL ENERGY DISTRIBUTION OF THE CRAB 
NEBULA 



The com pilation of observationa l data used here is sum- 
marized by lAharonian et alj (|2004|) and re ferences therein. 
The radio data ( Macfas-Perez et all 1201 fj) have been cor- 
rected for a secular decline of -0.18 %/yr to a common date 
(01/01/2000). Additionally, new data are added and listed in 
Table [I] while the entire compilation of data is displayed in 
Fig. |l(a)| and |(b)| The far infrared (FIR) observations from 
Spitzer, ISO, and Scuba deviate from a simple power-law ex- 
trapolation of the radio spectrum. The flux of the o ptical line 



emission is taken from [Davidson & 



le frux or the optical line 
Fesenl (1 19851) . bavidsonl 



(119871) and lHester et alj (Il990b . The optical line emission of 
the filaments in the nebula is estimated in the following way. 
The high-resolution spectral observations of individual fila- 
ments have been corrected for extinction dHesteretal.ll 19901) 
and scaled to match the global emission from the filaments 
(see the discussion in lDavidson & Fesenlll985l) . 

New X-ray observations by the XMM-Newton and INTE- 
GRAL observatories have been chosen to replace older mea- 
surements. Both, the XMM-Newton and INTEGRAL (with 
the instruments SPI and IBIS/ISGRI) observatories are cali- 
brated on th e basis of detailed simulations and labora tory mea- 
surements dAttie et al.ll200l iGondoin et al.lll998albl) . This 
approach differs from commonly used methods in which cor- 
rections of the instrument's response function are applied to 
reproduce a specific spectral shape and flux of the Crab Neb- 
ula. As a result, corrected measurements are model-dependent 
and therefore we chose not to include them here. The cameras 
of XMM-Newton spatially resolve the Crab Nebula, whereas 
the measurements of SPI and IBIS/ISGRI may include con- 
tributions from the pulsar, possibly leading to higher fluxes in 
comparison to the XMM-Newton observations. We stress that 
the difference in flux normalization between XMM-Newton 
and SPI are beyond the systematic errors quoted. To combine 
the two measurements, we chose to scale the flux of the SPI 
data by a factor of 0.78 to match the extrapolation of the mea- 
surement of XMM-Newton. 



A. Constant B-field model 

The nebula is assumed to be filled with relativistic electrons 
with an averaged total differential number dA^i/dy. To explain 
the change in the continuum between the radio and infrared 
(see above), we distinguish between radio and wind electrons. 
The two spectra are assumed to follow power laws with appro- 
priate cut-offs. For the radio-emitting electrons, a power law 
between yj. = 22 (a value only constrained towards higher 
values of y and kept constant in the fit), and y r mdS 
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is sufficient, while for the wind electrons a broken power law 
with a superexponential cut-off at the lower end is required: 
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As assumed by iHillas et all dl998l) . the electron-number 
volume density nZ(y, r) in the nebula is taken to drop off radi- 
ally following a Gaussian function such that 



n d(r. r) = exp 



2pHy) 



dy 



(3) 



The width of the Gaussian p(y) decreases with increasing 
Lorentz factor y to account for the observed shrinking of the 
nebula (see Appendix lAl for further details). The volume of 
the nebula is assumed to be filled with an entangled magnetic 
field of constant field strength. Within the nebula, various seed 
photon fields are upscattered by the electrons via the inverse 
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Compton process. The effective density of seed photons « see d 
is calculated by convolving the electron density with the pho- 
ton density and summing the contributions from several differ- 
ent components: (1) synchrotron radiation, (2) emission from 
thermal dust, (3) the cosmic microwave background radiation 
(CMB), and (4) optical line emission from the nebula's fila- 
ments. Similar to the electron population, the spatial photon 
densities are approximated by (energy-dependent) Gaussian 
distributions, whereas the photon density of the CMB is as- 
sumed to be constant throughout the nebula. The relevant re- 
lations for the seed photon fiel ds are summarized in Appendix 
lAl (for further details see also iHillas et al.|[l998h . The result- 
ing synchrotron and inverse Compton emission (at frequency 
v) is found by 



f*y™±(j* +J *). (4) 



The specific single particle emission functions X. v for syn- 
chrotron (Sy) and inverse Compton (IC) processes are given 
in AppendixlBl 

For a given value of the magnetic field, the shapes of the 
two electron spectra are varied until the resulting synchrotron 
spectrum reproduces the observational data. The best-fit val- 
ues for the ten free parameters describing the electron spec- 
trum at a fixed magnetic field strength are determined by 
means of a least-squares algorithm (Levenberg-Marquardt), 
which provides the closest match of the expected synchrotron 
spectrum (between 1 and 10 14 GHz) with the observational 
data (see Table HIl for the results of the fit and the discussion 
in Section 

[IIIJ- The resulting ^/d.o.f = 214.5/217 is be- 
low unity after a relative systematic uncertainty of 7 % was 
added in quadrature to the statistical error quoted for the data 
in the literature. The minimization procedure is not sensitive 
to the particular choice of starting values and it converges reli- 
ably. The resulting covariance matrix allows for an analysis of 
the correlations between various parameters. As expected, the 
normalizations are anti-correlated to the power law indices. 
Additionally, there are modest anti-correlations between the 
position of the break in the spectrum of the wind-electrons 
with the power-law indices. The matrix of correlation coeffi- 
cients is listed in AppendixICl 

The predicted inverse Compton component above 700 MeV 
is then used to calculate x 2 (B) f° r a range of B-field values 
using only the Fermi data. The best-fit value of B and its 
statistical uncertainties are estimated by adjusting a parabola 
to x 2 (B) an d calculating its second derivative. After taking 
the systematic energy uncertainty on the global energy scale 
of the Fermi/LAT data into account, AE/E = ^™ (see e.g. 
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lAbdo et alj|2009l) . the average Z?-field is found to be 
B = (l24 + 6 (stat.) + l b 5 (sys.)) yuG. 



(5) 



The reduced ^ 2 /d.o.f. = 6.37/13 ~ 0.49 for the fit of the 
inverse Compton component to the Fermi data indicates that 
the statistical errors on the Fermi differential flux may be 
slightly overestimated. 



The resulting broad band SED is displayed in Fig. |l(a)| ( 
solid and dashed blue curves) including synchrotron and IC 
emission and thermal emission from the dust in the nebula 
and optical line emission from the filaments (dashed magenta 
line). The inverse Compton component including the various 
contributions of different seed photon fields is shown in de- 
tail in Fig. [2] For convenience, an analytical parametrization 
of the derived energy flux above 1 GeV is presented in Ap- 
pendix [D] This parametrization is especially useful to derive 
cross calibration factors as introduced in Section[IV]for other 
or future instruments that measure the same energy range. 

For the thermal dust emission a graybody spectrum is used. 
By fitting the combined spectrum (thermal and nonthermal 
emission) to the data a temperature of T = 93 K is derived. 
The graybody pe aks at (vf v ) mj , x 
With the relation (IGehrz et alll998l) 



6.45 x 10 9 ergs s 1 1 



-2 



M dust 



1.36(v/ v ) max 

ap And 2 
3£r SB r 4 



(6) 



the dust mass Mdust can be estimated. Here <x SB denotes the 
St efan-Boltzmann con stant. We adopt the same assumptions 
as lTemim et al.l d2006l) : i.e., we assume a dust grain size of a = 
10 /jm and graphite grains with a density of p — 2.25 g cm ~ 3 . 
This leads to a dust mass of Mdust = 4x 10 ' 4 M H , about 40% 
of the value obtained bv lTemim et alj d2006l) . 



B. MHD flow model 

While the treatment presented above provides an accurate 
match of the measured broad band SED, a more physical de- 
scription of the evolution of the injected particles and mag- 
netic field in a spherical volume has been suggested by KC84. 
In addition to the synchrotron component already studied by 
KC84, AA96 extended this approach by investigating the 
high-energy (inverse Compton) component of the SED. The 
MHD solution requires the shock distance r s , the magnetiza- 
tion parameter cr, i.e. the ratio of the Poynting flux to the par- 
ticle energy flux at the position of the shock, and the shape of 
the injected electron spectrum (particle number per unit vol- 
ume n(y) in the interval y to y + dy) as main input parameters. 

Previously, RG74 assumed a shock distance of r s = 
0.10 pc, which in the meantime has become the canonical 
value, whe reas recent high spati al resolution Chandra ob- 
servations dWeisskopf et al.ll2000h indicate that the shock (if 
identified with the bright inner ring in the X-ray image at 14" 
distance to the pulsar) resides at a distance of r s = (0.14 ± 
0.01) pc. 

Based upon the initial analysis of KC84, the magnetiza- 
tion paramater was considered to be less than 1% with typ- 
ical values ranging from cr = 0.001 to cr = 0.005. With 
such low values of cr, both the observed break in the spec- 
trum between optical and radio wavelengths and the general 
morphology of the nebula are reproduced for the most part. 
However, the asymmetry in the brightness of the far and near 
sides o f the torus observ ed in X-rays and the measured flow 
speed dMori et al. | |2004 imply a higher value of cr between 
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(a) The SED of the Crab nebula calculated in the framework of the constant B-field model. The open blue data points have been included in the 
fit for the synchrotron part and the filled blue points used to determine the best-fitting magnetic field, 
r 1 . 1 . 1 1 . 1 . 1 . . 1 1 1 . . 1 1 1 . q 




FIG. 1. The SEDs for the constant B-field and MHD flow model. The data compilation is based upon Aharoni an et alj d2004l) with updates 
listed in Table U and modifications (e.g. scaling of the SPI flux) described in the text. See Section [TT] for further details. The synchrotron 
emission fitted to the data is shown by the blue solid line while the remaining components (low energy part and inverse Compton emission) 
naturally result from the model. Therefore, blue data points have been included in the fit whereas the gray points have been excluded. In 
addition to the nonthermal continuum emission, the contribution of dust and line emission from the filaments is indicated with a magenta 
dashed line. The lower panels show the residuals of the fit. . 

0.01 and 0. 13, more likely close to <x = 0.05 dMori et alj2004t geneously in the nebula with 
IShibata et alJl2003l) . 

Following the approach suggested by AA96, the broad band "radio(y) = my~ Sr exp I — J— J , (7) 

SED is calculated (see Eq.|H including the same external seed Tmax/ 
photon fields as described above and two populations of elec- 
trons. The radio-emitting relic electrons are distributed homo- whereas a population of wind electrons is injected at the 
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log 10 (E/eV) 

8 9 10 11 12 13 14 




23 24 25 26 27 28 

log 10 (v/Hz) 



FIG. 2. The Fermi/LAT and IACT data points are shown along with the total IC flux of the constant B-field (black line) and individual 
components from the different seed photon fields (constant B-field model): (1) synchrotron, (2) thermal dust, (3) CMB, and (4) line emission 
from filaments. 



shock: 

«wind(r) = qo(y + rl n )^" «p(— J- 1 . (8) 

The radiative and adiabatic cooling of the wind electrons is 
treated in the same way as suggested by AA96. The best- 
fitting value for the seven parameters describing the electron 
spectra are found in a similar way to what is described above. 
For a fixed value of cr, the parameters describing the elec- 
tron distribution (« , S r , y' m . dx , q , yj^, S w , y^ ax ) are varied, 
until the predicted synchrotron emission is matched best to 
the same data as used for the constant B-field model. The 
procedure is repeated for a range of values of cr until an 
absolute minimum is found at cr = 0.0045 ± 0.0003 and 
X 2 /d.o.f. = 397.02/222 * 1.79 with a probability of obtain- 
ing a higher value of x 1 per chance P(> x 2 ) — 4.8 x 10~ 12 (see 
TableHTTlfor the best-fit values and Appendix[C]for the correla- 
tion coefficients). The high^ 2 value indicates that the simple 
MHD-flow model fails to describe the synchrotron part of the 
SED in detail, while the overall shape is certainly correct. It is 
noteworthy that the high value of y' m . dx implies that most of the 
optical emission is produced by the same population of elec- 
trons as are responsible for the radio emission. In this case, 
the spatial extension predicted at optical frequencies would be 
similar to the extension of the radio nebula which clearly con- 
tradicts observations. When looking at the inverse Compton 
component calculated in this model, the agreement between 
the Fermi part of the spectrum and the model is fairly good, 
even though the shape of the synchrotron cut-off is harder than 
the measured spectrum. At energies beyond the position of the 



peak in the inverse Compton part of the SED, the model spec- 
tra are considerably harder than the actual measurements. The 
discrepancy is directly related to the mismatch that is evident 
in the X-ray part of the synchrotron spectrum. 



C. Comments on the two different model approaches 

The average magnetic field derived above (see Eq.|5J is dis- 
played as a yellow band in Fig. [3] together with the MHD so- 
lution of the magnetic field for cr = 0.0045. It is obvious that, 
in a spatially varying magnetic field, the effective magnetic 
field seen by the population of electrons radiating at different 
gamma-ray energies varies as well. Therefore, it should be 
possible to determine cr from the combined spectral measure- 
ment of the synchrotron and inverse Compton components. 
However, it appears that the broad band spectrum in the KC84 
model is not consistent with the data, which implies that some 
of the assumptions may have to be refined before using the 
model further in a quantative way in combination with the im- 
proved data available. 

A straightforward extension of the spherical KC84 model 
would be to incorporate an a symmetric flow with a mod- 
ulated magnetization (see e.g. iDel Zanna et all 120041 120061: 
iKomissarov & Lvubarskvll2004t IVolpi et al.ll2008l) . Depend- 
ing on the particular way the magnetization just upstream of 
the shock varies, the superposition of the emission from re- 
gions with a different downstream flow magnetization could 
be arranged to be closer to the observations than the single 
cr-model of KC84. In fact, by e.g. superposing the emission 
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of two separate regions with <x a 0.01 and cr * 0.001 would 
provide a better overall description of the data. Given that the 
resulting volume-averaged magnetic field in the downstream 
region could approach an effectively constant field, the result- 
ing inverse Compton emission would be comparable to the 
simple model considered here. 

III. ELECTRONS IN THE NEBULA 

The total number electron spectrum dA^i/dy of the constant 
B-field model is shown in Fig.|U 

A. Radio electrons 



Lorentz factor of the radio electrons as induced only by ra- 
diative cooling (ignoring adiabatic losses), an upper limit on 
the magnetic field can be derived: 

B< 385( A ) ,uG. (9) 

1 1.74 x 10 5 / \956yr/ p 

A lower limit on the magnetic field can be estimated by vary- 
ing the magnetic field and the normalization Nq of the radio 
electrons until the inverse Compton component overshoots the 
Fermi data. The constraint on the magnetic field is estimated 
to be B > 100 fiG. The two limits are consistent with the value 
derived for the constant B-field model as well as with the pre- 
diction by the KC84 model. 



lAtovanl d 19991) has suggested that the radio-emitting elec- 
trons (see Section UlAl and Fig. |4]i were injected in the phase 
of rapid spin-down during the initial stages of the pulsar- 
wind evolution. Observations of time-variable emissions from 
radio wisps with a hard spectrum i ndicate ongoing accel- 
eration of radio emitting electrons dBietenholz et al.l l2004t 
iBietenholz & Kronberg|[l990l: iKronberg et al.lll993l) . How- 
ever, the injection rate determined from the observations of 
the wisps is not high enough to explain the total population 
of radio-emitting electrons, therefore we consider them to be 
relic electrons. 

The radio electrons loose energy radiatively and adiabat- 
ically as the nebula expands. Interpreting the maximum 
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FIG. 3. The black horizontal line and the yellow shaded region mark 
the volume averaged magnetic field and its uncertainties as derived 
from the constant B-field model (see Eq. The magnetic field 
configuration as derived from the MHD solution of the downstream 
flow (KC84) is shown for <x = 0.0045 for a shock distance between 
0.13 pc < r s < 0.15 pc, where r s = 0.13 pc results in the upper and 
r s = 0. 15 pc in the lower boundary, respectively. The dashed-dotted 
lines mark the relative emissivity (normalized to z = 1) at specific 
energies. The horizontal line at 100 fiG indicates a lower limit on the 
B-field, see Sectionlllll 



B. Wind electrons 

The continuously injected wind electrons produce the bulk 
of the observed SED above sub-mm/FIR wavelengths via syn- 
chrotron emission. The spectrum is given in Eq. [2] The ra- 
diatively cooled spectrum of the wind electrons has a spec- 
tral index of S w = 3.23 = 2.23 + 1, which can be ex- 
plained naturally by 1 st -order Fermi acceleration at an ultra- 
rela tivistic shock with subsequent synchrotron cooling (see 
e.g. lKirk&Duffvlll999h . The low-energy cut-off of the wind 
electrons is found to be y w . = 4.24 x 10 5 , which is an- 

• mm 

ticipated in the model suggested by KC84. It requires the 
average Lorentz factor (y w ) of the isotropized (downstream) 
electrons to be similar to the bulk Lorentz fac tor y„ of the 
upstream electrons, (y w ) - 7.33 x 10 5 » y, ( Arons|[T996t 
iKennel & Coronitilll984tlKundt & Krotschecklll980l) . An ad- 
ditional feature present in the hard X-ray spectrum, which 
shows a softening at ~ 130keV, corresponds to a break with 
AS = 0.43 at y* . = 3.01 x 10 8 in the electron spectrum. The 
origin of this feature in the electron spectrum is very likely 
related to the injection/acceleration, given that it can hardly 
be related to energy-dependent escape. (The X-ray emitting 
electrons suffer cooling well before escaping the nebula.) The 
value of AS = 0.43 could hint at an energy dependent effect 
similar to diffusion in a Kolmogorov-type turbulence power 
spectrum. 

The total energy of the radio and wind electrons, respec- 
tively, is found to be 



, r div*. 

J r 17 



dy = 3.10 x 10 48 ergs, 



(10) 



dm 



E w = mc 2 j y^- & 
i 



dy 



indicating that the total energy in electrons is much less than 
the energy released through the spin-down of the pulsar. That 
both relic electrons and wind electrons have rough ly equal 



energy is consistent with the expectation given by lAtovanl 
(1999). 
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E[eV] 




FIG. 4. The two components of the electron spectrum used to calculate the broad band emission of the Crab Nebula in the constant B-field 
model. Red solid line: radio electrons; blue solid line: wind electrons. The black dashed lines indicate the values fot the minimum, maximum, 
and break energies. 



TABLE II. Parameters for the electron spectrum used for the constant B-field model. At the energy y rain the radio electron spectrum cuts off 
sharply whereas the wind electrons cut off superexponentially. The energy y max denotes a sharp cut-off for the radio and wind electrons. 







Constant B-field model 


Parameters 




Radio 


Wind 


Normalization constant 


ln(7V ) 


120.0(1) 


78.6(3) 


Low energy cut-off 


ln(7n,in) 


3.1 


12.96(3) 


Super exponential cut-off parameter 


/>' 




2.8(4) 


Break position 


hl(7break) 




19.5(1) 


High energy cut-off 


hlCymax) 


12.1(7) 


22.51(3) 


Spectral index 


s 


1.60(1) 


3.23(1) 


Spectral index (after break) 


S + AS 




3.76(3) 



TABLE III. Parameters for the electron spectrum for the MHD flow model. Compare Eqs.[7]and[8] 







MHD flow model 


Parameters 




Radio 


Wind 


Normalization constant 


In0i ), ln(go) 


-11.41(2) 


-1.32(1) 


Low energy cut-off 


ln(y rain ) 




13.94(2) 


High energy cut-off 


ln(y raax ) 


13.55(3) 


22.60(1) 


Spectral index 


S 


1.58(1) 


2.32(1) 


Magnetization parameter 


<T 


0.0045(3) 



IV. CROSS CALIBRATION OF IACTS & FERMI from the Crab Nebula in the high-energy regime. 



The updated model of the SED of the Crab Nebula provides 
an opportunity for a cross calibration between ground-based 
air shower experiments and the Fermi/LAT. The method is 
demonstrated here with the imaging air Cherenkov telescopes 
HEGRA, H.E.S.S., and MAGIC but is generally applicable 
to any other experiment that measures the flux and spectrum 



In general, the energy calibration of IACTs is done in- 
directly with the help of detailed simulations of air show- 
ers and the detector response. However, the remaining sys- 
tematic uncertainty on the absolute energy scale typically of 
15% leads to substantial differences in the observed flux and 
position of cut-offs in the energy spectra between different 
IACTs and also between Fermi/LAT and IACTs. First at- 
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tempts to cross-calibrate the I ACTs among each other have al- 
ready used the Crab Nebula dHorns et alj|2005l IZechlin et al.1 
120081) . Since the known gamma-ray spectra lack sufficient 
sharp features, an absolute cross calibration with lines, etc., is 
not feasible. So far, efforts to cross-calibrate Fermi/LAT with 
IACTs have focused on using the overlapping energy range 
dBastieri et al.ll2005l) . Cross calibration between Fermi/LAT 
and IACTs indirectly provides a means of benefiting from 
the careful beam-li ne calibration of the Fermi/LAT (see e.g. 
lAtwood et ai]l2009h . 

The cross calibra tion is accomplished in the following way 
dMever et alj|2009l) : The average magnetic field of the model 
is adapted to the Fermi observations. The statistical uncer- 
tainty on the B-field in Eq. [5] translates into statistical errors 
on this reference model (see Table llVI ). For each IACT an en- 
ergy scaling factor siact is introduced to correct the measured 
energy E meas to a common energy scale E such that 

E = E meas ■ siact- (12) 

The scaling factor siact for each instrument is determined via 
a ^-minimization in which the energy scale is changed ac- 
cording to the formula above until the data points reproduce 
the model best. The resulting scaling factors for the different 
instruments are listed in Table [IV] with the statistical errors 
and the reduced minimum^ 2 -values before and after the fit. 
To illustrate the result, Figs. |5(a)| and |5(b)| compare the un- 
sealed and the scaled data points with the model. It is evident 
from these figures that the scaled data points fit the model 
better. This is also quantified by the ^-values before and af- 
ter the application of the scaling factors. All scaling factors 
lie within the aforementioned 15% energy uncertainty of the 
IACTs. The complete SED with the scaled data points and the 
model calculations of SectionHDis shown in Appendix[E] 

The cross calibration eliminates the systematic uncer- 
tainty of the energy scale of the IACTs and adjusts them 
to a common one shared with the Fermi/LAT. However, the 
Fermi/LAT's absolute energy uncertainty remains, but it im- 
plies an improvement from +15% to +5% and -10%. 

As a first application of the cross calibration, we de- 
rive u pper limits on th e diffuse y-r ay background. Both 
Fermi dAbdo et al.l2009l) and H.E.S.S. dAharonian et al.l2008l 
120091) have measured the cosmic ray e + e + spectrum. Un- 
like Fermi/LAT, the telescopes from H.E.S.S. cannot accu- 
rately distinguish between showers induced by electrons (or 
positrons) or photons, such that up to » 50% of the observed 
electromagnetic air showers could be induced by photons. 
Therefore, H.E.S.S. actually measures electrons and diffuse 
background photons. Taking the difference of the two mea- 
surements we can derive an upper limit on the intensity of the 
y-ray background. The scaling factors derived above are now 
used to convert the H.E.S.S. data to the Fermi/LAT energy 
scale, which substantially reduces the systematic uncertainty 
on the observed intensity, given that the electron spectrum fol- 
lows a soft power law with oc E~ 3 . An important result of the 
cross c alibration is that th e peak in the spectrum observed by 
ATIC dChang et alj|2008l) appears more unlikely after apply- 
ing the scaling factors. 

The upper limits are derived by taking the difference of the 



two measurements from the energy region covered by both in- 
struments. This corresponds to the first six H.E.S.S. points of 
the low energy analysis in Fig. [6] The remaining systematic 
errors are taken into account for deriving the upper limits: the 
flux points of the H.E.S.S. measurements are shifted to their 
maximum value allowed by the systematic uncertainties while 
the Fermi points are shifted to the minimum value. This gives 
a conservative approximation for the upper limits. In general, 
the scaling factors can be used, e.g., to improve contraints on 
dark matter mo del parameters that rely on comb ined measure- 
ments (see e.g. lKistler & Siegal-Gaskinsll2010h . 



V. SUMMARY 

Updated models and data compilation for the SED of the 
Crab Nebula have been presented. The MHD flow model 
based upon KC84 and AA96 does not provide a satisfactory 
description of the available data and requires refinements of 
the underlying assumptions, e.g., relaxing spherical symmetry 
in axisymmetric numerical calculations as carried out by Volpi 
et al. (2008). A straightforward modification is, among oth- 
ers, the introduction of an anisotropic wind with variations of 
cr when moving out of the equatorial plane towards the polar 
regions of the outflow. The simplified approach of a constant 
magnetic field pursued here has the benefit of a smaller total 
number of parameters, even though the prescription of the in- 
jection spectrum in the MHD model is simpler (7 instead of 
10 parameters) and physically more meaningful. Most impor- 
tant for the task of cross-calibrating the instruments, the in- 
verse Compton component predicted in this model accurately 
describes all observational data above 1 GeV for a magnetic 
field of B - (124 ± 6 (stat.) ^ (sys.)) //G. The comprehen- 
si ve SED allow s for an estimate of the dust mass. In contrast 
to iTemim et alj d2006l) . we obtain a value about 40% lower, 
namely Md us t = 4 x 1O~ 4 M . The radio (relic) electrons pro- 
vide additional (independent) constraints on the average mag- 
netic field in the entire radio nebula. Using the endpoint of 
the electron spectrum the magnetic field is constrained to be 
smaller than « 385 yi/G, while the Fermi/LAT observations set 
a lower limit at approx. 100 //G. 

The model describes the broad band data at high en- 
ergies well enough to derive energy scaling factors for 
cross-calibrating the Fermi/LAT instrument with a variety of 
ground-based instruments. The cross-calibration eliminates 
the systematic uncertainty of the different energy scales used 
by the IACTs and ultimately reduces the global uncertainty to 
the Fermi/LAT calibration uncertainties of +5% and -10%. 
An application of the cross calibration were presented, and up- 
per limits on the diffuse photon background between 200 GeV 
and 1 TeV derived by combining Fermi/LAT with H.E.S.S. 
measurements. The excess measured by the ATIC collabora- 
tion seems unlikely with the scaled H.E.S.S. observations. 
The presented cross calibration can be used universally when 
combining observational data from various instruments, as 
well as in cases where the absolute energy scale is impor- 
tant. Specifically for soft spectra with cut-off features, such 
as in observations of objects at cosmological distances suffer- 
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(a) The IC model (constant B-iield model, solid black line) with (b) The same curve as in |5(a)|bu t with the scaling factors of Eq. [T2l and 

measurements from IACTs and Fermi/LAT. No energy scaling is Table llVl applied to the data, 

applied. 



FIG. 5. Comparison between the unsealed and scaled data of the IACTs. 



Instrument Scaling factor stact Stat, error As ^ cfore /d.o.f. l y^ fter /d.o.f. 



Fermi/LAT 


1 


+0.05 -0.03 




0.49 


HEGRA 


1.042 


±0.005 


7.652 


1.046 


H.E.S.S. 


0.961 


±0.004 


11.84 


6.476 


MAGIC 


1.03 


±0.01 


1.671 


0.656 



TABLE IV. Energy scaling factors of the IACTs for the cross calibration. 



H.E.S.S. - sys. error 




E [GeV] 



FIG. 6. e~ + e + spectrum reported by H.E.S.S. and Fermi. The cross 
calibration factor is applied to the H.E.S.S. energy scale, so the un- 
certainty on the global energy scale is reduced. The abbreviation 
l.e.a. stands for low-ene rgy analysis as applied to the H.E.S.S. data in 
lAharonian etail ll2009h. For comparis on, the PAMELA data points 
are taken from lMocchiutti et al. (2010). 



ing from absorption on the extragalactic background light, the 
cross calibration can be useful for providing more stringent 
upper limits. 



Appendix A: Spatial distribution of seed photons and electrons 
in the nebula 

To calculate the seed photon fields required to calculate the 
IC flux, it is neccessary to convert the observed photon flux 
to the corresponding photon densi ties. This is done by fol- 
lowing the approach suggested by iHillas et al.l d 1998b - Both 
the photon densities (apart from the CMB contribution) and 
electron density are assumed to follow Gaussian densities in 
the distance r from the Nebula's center, exp[-r 2 /(2cr 2 )] and 
exp[-r 2 /(2p 2 )], respectively. The variances a 2 and p 2 are 
energy dependent and estimated from ob servations. For the 
Gaussian describing the photon density IHillas et ail d 1 998b 
find 

arctanf — ) = 
\d) 

'3.93xl(T 4 fore < 0.02 eV, 

(Al) 

^0.47 + 3.46 [oof^v] ' 09 ) x 10 ~ 4 for e - °- 02 eV - 

Assuming that the photons are produced only via synchrotron 
radiation in a uniform magnetic field, the variance of the elec- 
tron density follows directly from Eq. IB3l with the aforemen- 
tioned averaging over 8, 



arctan 




10 



3.93 x 10" 



(o.47 + 28.9r- 017 [ T o^Gr°' 09 )xlO- 



fory < 6.65 x 10 4 

fory > 6.65 x 10 4 . 

(A2) 



Only those photons and electrons with overlapping distri- 
butions can interact, so that the total distribution is found by 
convolution to give a total variance of <x 2 + p 2 . For a photon 
production rate S v , which is the sum of the contributions from 
synchrotron radiation, thermal dust emission, and optical line 
emission, we find 



"Seed - 



4ttc(p 2 + cr 2 ) 



«CMB ■ 



(A3) 



The photon production rates are calculated using the individ- 
ual luminosities, S v = L v /h, where 



£ v Sy d T , 



sync _ f d^Vd 

Jv "J dy 
i 



(A4) 
(A5) 



The fluxes from the line emissions are approximated by 6- 
distributions, and the L, are given in the references, see Sec- 
tion|IIl The luminosity for the dust emission is given by a gray 
body spectrum and the photon density of the CMB is calcu- 
lated from a black body with a temperature of T — 2.726 K. 
The total seed photon field can then be used to calculate the 
inverse Compton emissions by means of Eqs. IB2l andl4l 



the IC distribution function fic can be written as 

/icO,v,y) = 

1 f4ey^/(mc 2 )l 

2q\nq + (1 + 2q)(l - q) + - \ \ Ji (1 - #5) 

2 1+ 4-eyq / (mc z ) 



Appendix C: Correlation between the fit parameters 

The covariance matrix has been calculated at the position 
of the minimum in the ^-function. The diagonal elements 
of the covariance matrix have been already listed in the form 
of estimates of the error of the individual parameters, see Ta- 
bles [II] and [III] Here, we list the correlation matrix defined 
as cor(i, j) = -y/cov(/, j)/(o-jcrj) in Table [V] for the constant 
B-field model and in Table EDfor the MHD flow model. 



Appendix D: Parametrization of IC flux 

For convenience, the energy flux vf v of the IC emission for 
1 GeV < E < 1 PeV is parametrized as a function of energy, 
namely a 5 th -order p olynomial in double-l ogarithmic repre- 
sentation (as done bv lAharonian et al.ll2004l) : 



logi 



Vfv 



ergs 



5 , 
= Yj Pi 1Og 10 ( 



— )• 

TeV/ 



(Dl) 



Appendix B: Single-particle emission functions 

The single-particle emission functions for particles of 
mass m and energy ymc 2 for synchrotron _£y y and in- 
verse Compton emis s ion X ?, c used in Eq. [4] are taken from 
iBlumenthal & Gouldl(ll970b : 



----- — 



V3e 3 Bsin# 



CO 

-f 



v/v c 



hv 



4 c = ^»v j d 



"seed(e) , , . 

e /ic(e, v,y). 



(Bl) 



(B2) 



hv/(4y 2 ) 



We denote the photon energy after scattering by hv and be- 
fore scattering by e, the Thomson cross section by <x x and the 
electron charge by e. The critical frequency v c is defined as 



V c = 



3e 
Amnc 



BsinOy 2 , 



(B3) 



and Ks/3 (x) stands for the modified Bessel function of frac- 
tional order 5/3. The electron pitch angle 8 is averaged where 
the value (sin 6} = V2/3 is adopted. Introducing the kine- 
matic variable q, 

hv 

q ~ 4ey 2 [l-hv/(ymc 2 )Y (B4) 



The coefficients are given in Table IVHI The relative error of 
the parametrization for the parameters p^, and ps are approx- 
imately 6%, less than 1% for po and p2 and about 1% for 
p\. The value of p4 is set to zero since its relative error is 
otherwise around 150%, and thus p4 is not neccessary for a 
satisfactory fit. 



Appendix E: Final SED 

Figure [7] summarizes the best fits for the constant B -field 
model and the MHD flow model (see Section [II]), together 
with all data points o f the references in Table U and in 
lAharonian et alj d2004l) . Likewise, the scaling factors for the 
IACTs introduced in Section|IV]are also applied. 
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TABLE V. The correlation coefficients between the fit-parameters for the constant B-field model. Since the matrix is symmetric, the lower 
trigonal part is not listed. 



c 

o r 


ln(N' n ) ln( 7 ,' mlx ) ln( 7 - n ) InOCJ ln( r ^) 


p 


a w 


O w T ZAO in^iVQ ) 


S r 1 


-1.00 -0.35 -0.02 


-0.03 


-0.02 


0.10 


0.05 


0.02 0.03 


HN r ) 


1 0.34 0.02 


0.03 


0.02 


-0.10 


-0.05 


-0.01 -0.03 


ln (rLx) 


1 0.24 


-0.05 


-0.03 


0.26 


0.08 


0.02 0.05 


ln(y w ) 

v/ mm 7 


1 


0.56 


0.13 


-0.57 


-0.80 


-0.25 -0.57 


ln (7break) 




1 


0.39 


-0.46 


-0.82 


-0.70 -1.00 


m (rmax) 






1 


-0.11 


-0.22 


-0.75 -0.38 


P 








1 


0.71 


0.19 0.47 


s w 










1 


0.42 0.84 


S W + AS 












1 0.69 


WW 












1 



TABLE VI. The correlation coefficients between the fit-parameters for the MHD flow model. Since the matrix is symmetric, the lower trigonal 
part is not listed. 





\n(q ) lnO ) ln( r ; ilax ) ln^) ln(y;; iax ) 


s w 


Sr 


hi((?o) 


1 -0.02 0.23 0.79 


0.94 


-0.00 


0.00 


ln(no) 


1 -0.50 -0.03 


-0.02 


-0.01 


0.01 


m (7max) 


1 0.31 


0.19 


0.01 


-0.01 




1 


0.65 


0.01 


-0.01 


in(rLx) 




1 


-0.11 


0.11 


s w 






1 


-1.00 


s r 








1 



TABLE VII. Parametrization of the IC flux. The coefficients corre- 
spond to Eq. ID 1 1 
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